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Abstract 

We review the properties of chaoticity and coherence in Bose-Einstein con¬ 
densation and correlations, for a dense boson system in its mean-field repre¬ 
sented approximately by a harmonic oscillator potential. The order parameter 
and the nature of the phase transition from the chaotic to the condensate 
states are studied for different fixed numbers of bosons. The two-particle 
correlation function in momentum space is calculated to investigate how the 
Bose-Einstein correlation depends on the degree of condensation and other 
momentum variables. We generalize the Bose-Einstein correlation analysis to 
three-particle correlations to show its dependence on the degree of conden¬ 
sation. 
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1 Introduction 


As is well known, a fundamental assumption for the occurrence of Bose-Einstein 
correlation (BEC) is the presence of a chaotic source of identical bosons jl][2]- The 
Bose-Einstein correlation occurs in a chaotic source but not in a coherent source 

OllEllellzllS]. 

The properties of chaoticity and coherence are complementary attributes. Both 
chaoticity and coherence should be examined on equal footings in a single theoretical 
framework with the description of both the BE condensation and BE correlations. 
In such a unified framework, it is then possible to investigate not only the states of 
chaoticity and coherence, but also the transition from a chaotic state to a coherent 
state. How can the degrees of chaoticity or coherence be quantified? Is the transi¬ 
tion from a chaotic state to a coherent state a first-order with a sudden onset, or is 
it a gradual transition that is closer to a second-order? What is the relevant order 
parameter that best describes the transition? How does Bose-Einstein condensation 
quantitatively affects the two-particle and three-particle Bose-Einstein correlations? 

Questions of Bose-Einstein correlations and condensation arise not only in atomic 
physics [5][Tni|TI] but also in high-energy heavy-ion collisions [5]0|^ where pions 
are the most copiously produced particles. The use of two-pion Bose-Einstein cor¬ 
relations to probe the source coherence was proposed at the end of 1970s [12] 
The introduction of the “chaoticity” parameter A of BEC in pions is only a tool to 
represent experimental data. However, the experimental measurement of A is beset 
by the presence of many other effects such as particle misidentification, long-live 
resonance decay, final state Coulomb interaction, non-Gaussian source distribution, 
etc. mm- The explanation of the experimental A results remains an open question. 
In 1993, S. Pratt proposed a pion laser model in high energy collisions and studied 
the influence of pion laser on two-pion Bose-Einstein correlation function and the 
chaoticity parameter [13]. In 1998, T. Csdrgo and J. Zimanyi investigated the effect 
of Bose-Einstein condensation on two-pion Bose-Einstein correlations M- They uti¬ 
lized Gaussian formulas describing the space and momentum distributions of a static 
non-relativistic boson system, and investigated the influence of the condensation on 
pion multiplicity distribution. In 2007, G. Y. Wong and W. N. Zhang studied how A 
in Bose-Einstein correlations depends on the degree of Bose-Einstein condensation 
or chaoticity, for static non-relativistic and relativistic boson gases within a spherical 
mean-field harmonic oscillator potential jl5| . The model can be analytically solved 
in the non-relativistic case and be used in atomic physics [9|[T0|[n|. The limiting 
conditions and circumstances under which the parameter A can be approximately 
related to the degrees of chaoticity were clarified |15j . A similar study for cylindri¬ 
cal static boson gas sources was completed [16] and the chaoticity parameter A in 
two-pion Bose-Einstein correlations in an expanding boson gas model was recently 
examined m- The investigation of chaoticity and coherence was also carried out 
using a model of q-deformed oscillator algebraic commutative relations [l8] and the 
model of partial indistinguishability and coherence of closely located emitters Ca¬ 
in another related topic, initial conditions such as the color-glass condensate (GGG) 
with the coherent production of partons |20| in heavy-ion collisions may also lead 
to condensate formation m- 
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Recently, experimental investigation of the source coherence in Pb-Pb collisions 
at y/SNN = 2.76 TeV at the Large Hadron Collider (LHC) was carried out by the 
ALICE collaboration |22j . A substantial degree of source coherence was measured 
|22) using a new three-pion Bose-Einstein correlations technique with an improve¬ 
ment over past efforts [^. Earlier work on three-particle correlations 

were carried outin|6l[l3]l27l|28ll23|30l|3T][32l|33l|M]- 

A proper theoretical framework to study the above topics is the theory of the 
Bose-Einstein condensation and correlations in their own mean field potential (15) . 
We would like to review the essential elements here and examine further the related 
question of three-body correlations. 

2 Bose-Einstein Condensation for attractively In¬ 
teracting Bosons 

We seek a description of chaoticity in Bose-Einstein correlations through the consid¬ 
eration of Bose-Einstein condensation. Why is Bose-Einstein condensation relevant 
to Bose-Einstein correlations (BEC)? Glauber in many private communications and 
in his talk in QM2005 suggested that the consistent experimental observations of 
A < 1 may be due partly to the coherence of the pions in Bose Einstein correlations 
|35] . Furthermore, there have been major advances in Bose-Einstein condensation in 
atomic physics UllOlin]. In particular, the works of Politzer |9j, and Naraschewski 
&L Glauber [10] reveals that BE condensation and the BE correlations are intimately 
related. 

We envisage the possibility of the occurrence of a Bose-Einstein condensation 
in dense boson media of identical bosons with the following reasoning |15l ITFl ITT] 

1. Identical bosons with mutual attractive interaction generate a mean field po¬ 
tential, which depends on the boson density p(r) as [36] 

V(r) = -47r/(0)p(r) - ifiw , (1) 

where /(O) is the forward scattering amplitude, and a is the length scale that 
defines the spatial region of boson occupation. 

2. Therefore, for a given length scale a, the huj of the underlying mean-field 
potential increases with increasing density p of the produced bosons. 

3. The order parameter that determines the degree of BE coherence or chaoticity 
is T/hui. Thus the order parameter T/Tbuj decreases with increasing boson 
density. 

4. For a given temperature T at freeze out, a high density of produced bosons will 
lead to a lower value of the order parameter T/huj, which in turn will lead to a 
greater condensate fraction /o=Vo/V, where N is the total number of bosons 
and Nq is the number of bosons in the lowest state. A greater condensate 
fraction /o brings about a greater coherence in Bose-Einstein correlations and 
a reduction in the degree of chaoticity. 
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In high energy heavy-ion collisions when bosons (gluons or pions) are copiously 
produced within a small region in a short time interval, the density of the bosons 
increases as the collision energy increases. Following the above reasoning, general¬ 
ized to systems with differential transverse and longitudinal spatial distributions, we 
expect the occurrence of boson condensation in high energy heavy-ion collisions at 
some high collision energies. It is useful to examine the Bose-Einstein condensation 
for bosons in an exactly solvable model. 


3 Bose-Einstein Condensation for Bosons in a 
Spherical Harmonic Oscillator Potential 


We consider a system of N bosons in a spherical harmonic oscillator potential, 
which arises either externally or from the bosons’ own mean-fields. We study how 
the occupation numbers of different states change as a function of the temperature 
T, in relation to the oscillator frequency Tiui. Bose-Einstein condensation occurs 
when the occupation number Nq for the lowest state (the condensate state) is a 
substantial fraction of the total particle number N. The degree of coherence or 
chaoticity is quantified by the condensate fraction /o = Nq/N, which varies as a 
function of the order parameter T/fia;. 

In such a study, it is important to use the proper statistical ensemble |9|. In 
a grand canonical ensemble, we fix the chemical potential and the temperature 
T, and we allow the number of particles in the n-th single-particle state to 
vary. We obtain the average occupation number for the single-particle state n to 
be Nn = The square fluctuation of Nn is then given by 

((a+a„ - (a+o„))^) « iV„(iV„-F 1). (2) 


As the fluctuation of Nn in a grand canonical ensemble is of the same order as 
the occupation number itself, we cannot treat the lowest n = 0 state in the grand 
canonical ensemble. The lowest n = 0 state needs to be treated in the canonical 
ensemble with a fixed total number of bosons. 

It was shown however that while the lowest n = 0 state needs to be treated 
in the canonical ensemble, the n > 0 state can be treated in the grand canonical 
ensemble without incurring large errors l^. We shall follow such a description for 
the ensemble of N identical bosons in a spherical harmonic oscillator potential. In 
such a canonical ensemble for the lowest n—0 state but a grand canonical ensemble 
for the n>0 states, the total number of bosons is fixed and yields the condensate 
number condition 


iV = iVo+ 


Z 


1-z 


E 


g„ze 

1 — ’ 


(3) 


where z = is the fugacity of the system, g„ is the degeneracy number 

9n=('n+i)(n+2)/2 for the n-th single-particle level, and e„ is the single-particle 
energy in the spherical harmonic oscillator potential 


£n = {n + 2>/2)‘hoj. 


(4) 
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For a given N, equation Q contains only a single unknown, z, which can be solved 
as a function of the order parameter T/fiuj. The solutions of z for N=25, 500, 
1000 and 2000 are given in Fig. 1, and the corresponding condensate fractions 
/o = Nq/N are shown in Fig. 2. 



T/ ho) 

Figure 1: (Color online) (a) The fugacity parameter z satisfying the condensate 
number condition Eq. ([^ for different boson numbers in a spherical har¬ 
monic oscillator potential, as a function of the order parameter T /huj and (b) 
an expanded view in the z ~ 1 region. 

We observe in Fig. 1 that the fugacity parameter z is close to unity in the 
strongly coherent region at low temperatures. In fact, the fugacity parameter z 
at T—0 assumes the value z{T—0) = N/{N + 1). For a given boson number N, 
the fugacity z decreases very slowly in the form of a plateau, as the temperature 
increases from T = 0. The plateau region persists until the condensate temperature 
Tc is reached, and z then decreases rapidly thereafter. The greater the number of 
bosons N, the greater is the plateau region, as shown in Fig. 1(6). For example, 
for N — 2000 the value of z is close to unity for 0 < T/Tiuj < 11 in the plateau,. 

We note in Fig. 2 that for a given value of the total number of bosons N in 
the spherical harmonic oscillator potential, the condensate fraction /q is close to 
unity when the order parameter T/Uuj is below a limit, and this limit depends on 
N. We can plot the condensate fraction /o as a function of the order parameter 
Tjhu!. The functional form of fo{T) can be approximated by 

c = / 1 - [iT/nij)/iTjnu;)]^ for (T/Aw) < (Tjhu;), 

> \Cl(l/fV)^0 for (T/Aw) > (Tc/M- ^ 
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Figure 2: (Color online) Solid curves represent the condensate fractions fo{T), 
calculated with the condensate number condition Eq. ([^, as a function of T/hoj 
for different boson numbers in a spherical harmonic oscillator potential. The 
abscissa labels for the corresponding chaotic fraction fj'{T) = [l — /o(T)] are 
indicated on the right. The dashed curves are the fits to the solid curve results 
of /o(T) with the function 1— [(T/?iw)/(Tc/?iw)]^ of Eq. ([^ where the values of 
Tc/hco for different N values are listed in Table I. 


The results from the above one-parameter fit to fo{T) are shown as the dashed 
curves in Fig. 2, to be compared with the fo{T) calculated with the condensate 
configuration condition Eq. © shown as the solid curves. The values of Tc/fno 
that give the best fit to fo{T) for different N values are listed in Table I. 

The above results provide a comprehensive description for the transition from 
a chaotic state to a coherent state. Fig. 2 indicates that the transition from the 
completely chaotic state with fo—0 to the state of coherence with /q—>-1 is a second- 
order-type transition under a gradual decrease of the order parameter T/huj. It is 
not a first-order phase transition. 

Table I. Critical order parameter Tc/Tnv obtained from (i) fitting /g as a function 
of T/hoj with Eq. (|^, and from (ii) the analytical formula of Eq. (|^. 


Number of Bosons 
N 

Tc/fiio obtained 
from fitting /g with Eq. (|^ 

Tc/hu) obtained 
with Eq.(|^ 

2000 

10.97 

11.00 

1000 

8.56 

8.53 

500 

6.63 

6.62 

250 

5.12 

5.13 


It is remarkable that the critical order parameter Tc/huj and the boson number 
N obeys the following simple relationship 

Tc/hu} = 0.6777iV°-3®®®®, (6) 


6 














as shown by the third column in Table I. Thus, the knowledge of N sufFices to 
determine the critical order parameter Tc/Tioj by the above simple equation and the 
knowledge of Tcjfiuj subsequently yields the approximate condensate fraction at all 
other temperatures by Eq. 


4 Single-particle and Two-Particle Density 
Matrices in Momentum Space 

The determination of the fugacity z from the condensate number condition 
allows the calculation of various physical quantities. Specifically, the one-body 
density matrix in momentum space is given by 

OO 

^'^Unip[)un{Pi){aian), (7) 

n—0 

where M„(p) is the single-particle wave function and the occupation number 
Nn={a\^an) can be inferred from the terms in the summation in Eq. ®. The 
two-particle density matrix in momentum space 

G^‘^HPuP2',P[,P2) = “k(PiX(P2)Wm(P2)Wn(Pl)(4«iaman) (8) 

klmn 

can be written in terms of one-body density matrices as mils] 

G'^'^\pi,P2;p'i,P2) = G^'^\pi,p[)G^'^Hp2,P2) + 

OO 

+ ^niPlhniP2)Un{P2)UniPl) 

n—0 




Figure 3: Two-particle distribution function expanded in terms of products of 
one-particle distribution functions in uncorrelated mean-field approximation. 

The uncorrelated part in the first two terms of the above two-particle density ma¬ 
trix, G^^'^{pi,p[)G''^\p 2 ,P 2 ) G^^\pi,P 2 )G^^Hp 2 jPi)< is represented schemati¬ 

cally by the diagram in Fig. 3. Our task is to obtain the correlated part arising from 
Bose-Einstein condensation represented by the last term in Eq. @. 

In the limit of a large number of bosons in a grand canonical ensemble for 
the non-condensed states, the contributions from the set of {n > 0} states in the 
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summation in Eq. can be neglected. We are left with only the n = 0 condensate 
state contribution for this summation. 

To describe the contribution from the n = 0 condensate state, we shall follow 
Ref. |5] [lU] and use the canonical ensemble which gives the canonical fluctuation 

((ajoo - (ajao))^) = (ajajaooo) - (aJao)(aoao) = 0(No). (10) 

Thus, we have 

(ajajooao) - 2(aJ)ao) (ajao) = - (dido) (dido) + 0(No). (11) 

In the limit of a large number of particles, we can neglect the last term 0{Nq) in 
the above equation which is small in comparison with the first term of order 
The two-particle distribution of Eq. is therefore 

G^'^Hpi,P2'^Pi,P2) = + |G(i)(p^,p2)|2 

-No\uo{Pi)\^\uo{P2)\^, ( 12 ) 

which gives the conditional probability for the occurrence of a pion of momentum 
Pi in coincidence with another identical pion of momentum P 2 - 


5 Two-Particle Momentum Correlation Function 

In BE correlation measurements, we normalize the probability relative to the proba¬ 
bility of detecting particle Pi and P 2 , and define the momentum correlation function 

G{Pi,P2) as 


GiPl,P2) = 


G^'^\Pi,P2-,Pi,P2) 


G(i)(Pi,Pi)G(i)(p2,P2)’ 

It is convenient to introduce the average and the relative momenta of the pair 


(13) 


P = (Pi +P2)/2, 


q=Pi -P2- 


(14) 


The momentum correlation function can be expressed as a function of the kinematic 


variables Pi and p^ or alternatively of p and q. From Eq. (121, we have the general 
expression for the correlation function 


G{p,q) = C{pi,P 2 ) = 1 + 


\G^^\Pi^P2)\^ - N^\MPi)\^\MP2)\^ 

G^'^HPi,Pi)G^^'>{P2,P2) 


(15) 


This is the general Bose-Einstein correlation function for all situations: coherent, 
chaotic, and the transition between coherent and chaotic systems. 


The evaluation of the correlation function C{p, q) in Eq. (151 requires the knowl¬ 
edge of G^^\pi^P 2 ) and the ground state wave function uq{pi). For a system of 
bosons in a spherical harmonic oscillator, the wave functions are all known, and the 
correlation function can be written out analytically. Specifically, we have 






( 16 ) 


Go{Pl,P2'-:kPhu}), 




Go{Pi,P2',t)= 


3/2 


^7r?i^(l — e 


exp 


(Pi+Pi)(coshr-l) + (pi-p2)' 

2 sinh r 


and the ground state wave function is 

/ a2 ^ 




.(17) 


(18) 


The knowledge of the single-particle {pi,P 2 ) and uo{p) will then allow the 
determination of the two-particle correlation function C{p,q). 

The correlation function C{p,q) in Eq. ( [T^ possesses the proper coherent and 
chaotic limits. For a nearly completely coherent source with almost all particles 
populating the ground condensate state, Nq N, the two terms in the numerator 
cancel each other and we have C{p, q) = 1, with the absence of the BE correlation. 
For the other extreme of a completely chaotic source with Nq-^N, the second term 
in the numerator proportional to Nq in Eq. (15| gives negligible contribution and 
can be neglected. The correlation function C{p,q) then becomes the usual BE 
correlation for a completely chaotic source. 


6 Evaluation of the Two-Particle Momentum Cor¬ 
relation Function 

For a given number of bosons iV in a spherical harmonic oscillator, the solution 
of fugacity z obtained as a function of the order parameter T/hui allows us to 
evaluate the momentum correlation function C(p,q) with Eqs. In Fig. 

4, we show C{p,q) for example for the case of iV=2000 for which the critical 
order parameter is Tc/huj—10.97 , as tabulated in Table I. We observe that the 
correlation function is a complicated function of the average pair momentum p 
and the order parameter T/huj. For p=Ti/a in Fig. 4(a), the correlation function 
C{p,q) at q—0 is close to unity for temperatures below and up to T/huj—9 (below 
Tc/huj), but increases to 2 rather abruptly at T/Tiuj=12, (above Tc/Tiui). For 
p=2fi/a in Fig. 3(5), the correlation function C{p,q) at q = 0 is substantially 
above unity and increases gradually as temperature increases. Forp=37i/a in Fig. 
4(c), the correlation function C{p, q) at q—0 is about 2 for all cases of temperatures 
examined. If one follows the standard phenomenological analysis and introduces the 
“chaoticity” parameter A to represent the intercept of the correlation function at 
zero relative momentum, then this parameter A is a function of the average pair 
momentum p and temperature T 

X{p,T) = [Cip,q = 0-,T)-l]. (19) 
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Figure 4: (Color online) The correlation function C{p, q) at different values of 
the pair average momentum pa/fi and temperatures. Figures (a), (b), and (c) 
are for p=l, 2, and 37i/a, respectively. 


We display explicitly the dependence X{p,T) as a function of p in Fig. 5(a) for 
different order parameters T/huj, for the case of = 2000. At T(fiuj=12, which is 
above the critical condensate order parameter of Tc/huj—10.97, the A parameter is 
1 for all p values. At T/huj—9, as p increases the A parameter rises gradually from 
^0.1 at p = fi/a and reaches the constant value of 1 at p—2.Ah/a. At T ITioj—6 and 
3, for which the systems are significantly coherent with large condensate fractions, 
the A parameter starts close to zero at p—fi/a, but as p increases the A parameter 
increases gradually to unity at p=2.9 and Z.lTi/a for TIfiuj—O and 3 respectively. 
The location where the A parameter attains unity changes with temperature. The 
lower the temperature, the greater is the value of p at which the A parameter attains 
unity. 

We conclude from our results that the parameter A(p, T) is a sensitive function 
of both p and T, and A(p, T) = 1 is not a consistent measure of the absence of the 
condensate fraction, as it attains the value of unity in some kinematic regions for 
significantly coherent systems with large condensate fractions at temperatures much 
below Tc. Only for the region of small p will the parameter A(p, T) be correlated 
with the chaotic fraction /t(T) of the system. 

It is interesting to note that experimentally measured values of A from different 
collaborations and different method of analysis [37l [38] exhibit an increase of A as 
Pt of the average momentum of the pair increase as shown in Fig. 5(b). There 
is a similar trend of increasing A as a function of pt- This may be an indication 
of the dependence of the correlation function on the average momentum of the 


10 


















Figure 5: (Color Online) (a) The parameter A as a function of p for different 
temperatures for A^=2000. (b) Experimental measured values of A as a func¬ 
tion of pt for AuAu Collisions at RHIC at ^snn= 200 GeV from the PHENIX 
Collaboration |37j and the STAR Collaboration |38) . 


pair arising for a partially coherent pion source. The increase of A as a function 
of the average pair momentum has also been obtained in the q-deformed harmonic 
oscillator model of Bose-Einstein correlations [18j . 


7 Bose-Einstein Condensation of Pions in their 
Men Fields 

With regard to heavy-ion collisions at RHIC LHC, it is instructive to raise the 
following question. If we have a pion system that has a root-mean squared radius 
ri-ms^lO fm, the number of identical pions N from a few hundred to a few thou¬ 
sand, at a freezeout temperature T—80 to 160 MeV, typical of those revealed by 
Bose-Einstein correlation measurements |371138] , then, what will be the condensate 
fraction /g? To answer this question, it is useful to calculate the root-mean-squared 
radius of the pion system as a function of the order parameter T/huj for a pion 
system with N—250 to 2000 as shown in Fig. 6. We can schematically represent 
the functional relation between rmis/a and T/huj in Fig. 6 as 

?’rms/a = fNir/huj). (20) 

For a given value of N and as a is equal to fij^/rrh^Ku:, the above equation 
contains only a single variable Tilo that can be determined as a function of T. 
Subsequently, the order parameter T/Tiuj and the condensate fraction /g can also 
be determined as a function of T as shown in Fig. 7. 

One finds that for the pion system with a given root-mean-squared radius of 10 
fm, the value of fiuj ranges from about 12 to 20 MeV for A^=2000 and about 20 to 
30 MeV for A^=250. The ratio of T/huj about 7 for Af=2000, and is about 4.5 for 
A^=250, as shown in Fig. 7(6). From these ratios of T/huj, one can use Fig. 2 to 
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Figure 6: (Color online) The root-mean-squared radius in unit of a and the root- 
mean-squared momentum in units of ?i/o, as a function of T/huj for different 
numbers of identical bosons in the system. 

find out the condensate fraction. The condensate fractions fo{T) for a pion gas at 
various temperatures with iV=2000 and N—250 are shown in Fig. 7(c). One finds 
that fo{T) is about 0.67— 0.8 for iV=2000 and is about 0.9 for N—250. 

We reach the conclusion from the above study that if a non relativistic pion 
system maintains a static equilibrium within its mean field, and if it contains a 
root-mean-squared radius, a pion number, and a temperature typical of those in 
high-energy heavy-ion collisions, then it will contain a large fraction of the Bose- 
Einstein pion condensate. For a relativistic pion system, while the absolute scale of 
the order parameter T/Uuj may change, the condensate fraction /o remains substan¬ 
tial |15| . Pion condensation will affect the parameter A in momentum correlation 
measurements. 


8 Three-particle Correlations and Coherence 

Bose-Einstein condensation has important influence on the three-particle correla¬ 
tion function. We can determine the dependence of the three-particle correlation 
function on the degree of Bose-Einstein coherence in a way similar to what has been 
carried out for two-particle correlations. 

The extraction of the coherence properties from experimental three-particle cor¬ 
relation data has the advantage that the problems of the resonances can be mini¬ 
mized. It has however the disadvantage that the statistics in the number of three- 
particle events may be lowered because of the restriction on the occurrence of 
three-particle coincidences. 

Recently there has much interest in three-particle correlation measurements 
[22] . Bose-Einstein condensation of pions in a heavy-ion collision may suppress 
Bose-Einstein correlations. Furthermore, initial conditions such as the color-glass 
condensate (CGC) with the coherent production of partons [20j may also lead to 
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Figure 7: (Color online) (a) the potential strength hoj, (b) the ratio T/huj, (c) 
the condensate fraction /q, and (d) the oscillator length parameter a for non- 
relativistic boson systems with N = 2000 and N = 250 in a static equilibrium 
with a Truis = 10 fm, plotted as a function of the temperature T. 


condensate formation m- Experimental results indicate the presence of a substan¬ 
tial condensate fraction |22) . It is of interest to formulate an analytical model to 
investigate how the three-particle correlation function will depend on the coherence 
of the underlying boson system. In a completely chaotic source when multi-particle 



Figure 8: Three-particle distribution function expanded in terms of one-particle 
distribution function in uncorrelated mean-field approximation. 

Bose-Einstein-type correlations are neglected, the three-particle correlation function 
can be written in terms of products of one-body distribution functions: 

G(3)(pi 

,P2-P3;p'i,P2,P3) 

= G(^)(1,1')G(i)(2,2')G(i)(3,3') + G(^)(1,2')G(i)(2,T)G(i)(3,3') 
-hG(i)(l,3')G(i)(2,2')G(i)(3,l') + G(i)(1,1')G(i)(2,3')G(i)(3,2') 
-hG(^)(l,3')G(^)(2,l')G(^)(3,2') + G(^)(1,2')G(^)(2,3')G(^)(3,1'), (21) 
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as represented by the diagrams in Fig. 8. With Bose-Einstein correlations, we can 
generalize our two-particle correlation case to the three-particle correlation functions 
and write down the three-particle correlation function as 


C{P1,P2,P3) 


G(3)(1,2,3;1',2',3') 

G(i)(1,10G(i)(2,2')G(i)(3,3') 




G(i)(l,2)G(i)(2,l)-jVoug(pi)zr§(p2) 

G(i)(l,l)G(i)(2,2) 

1) - N^ul{pi)ul{p3) 

G(i)(l,l)G(i)(3,3) 

G(i)( 2,3)GW(3,2) -fVgug(p2)ug(p3) 

G(i)(2,2)G(i)(3,3) 

GW(l,3)G(i)(2,l)G(i)(3,2)-iV3ug(pi)zr§(p2K(P3) 

G(i)(l,l)G(i)(2,2)G(i)(3,3) 

GW(l,2)GW(2,3)G(i)(3,l)-iVgug(pi)zrg(p2K(p3) . . 

G(i)(l,l)G(i)(2,2)G(i)(3,3) ’ 


The above correlation function C{pi,P 2 ,P 3 ) possesses the proper coherent and 
chaotic limits. For a nearly completely coherent source with almost all particles 
populating the ground condensate state, Nq^N, the terms in the numerator can¬ 
cel each other and we have C{pi,P 2 ,P 3 )=i, and the BE correlation is absent. For 
the other extreme of a completely chaotic source with No<tiN, the second terms 
in the numerators proportional to Nq give negligible contribution and can be ne¬ 
glected. The correlation function C{pi,P 2 ,P 3 ) becomes the usual BE correlation 
for a completely chaotic source. These results will allow the evaluation of the three- 
particle correlation function using the functions of {pi,P 2 ) and Uq{p) in Eqs. 
@-@. Different ways of re-combining some of the terms in Eq. ( [2^ in terms of 
two-particle correlation functions may allow one to extract quantities that minimize 
the systematic errors in two-particle correlation function measurements. 


9 Conclusions and Summary 

A proper framework to study Bose-Einstein correlations is the theory of Bose- 
Einstein condensation. We examine the condition for the occurrence of the Bose- 
Einstein condensation in an exactly solvable model. We place identical bosons in a 
spherical harmonic oscillator potential that arises either externally or approximately 
from its own mean fields. The order parameter is T/huj, the ratio of the tem¬ 
perature to the energy gap between the lowest and the first excited single-particle 
state. The degree of chaoticity or condensation is quantified by the condensate 
fraction /o = Nq/N which specifies the transition from a chaotic state to a coher¬ 
ent condensate state. The condensate fraction /o is a cubic function of the order 
parameter T/Tiuj. The critical order parameter Tc/Tiuj varies with the boson number 
N as Tc/T ujj= 0.6777 The transition from the completely chaotic state with 
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/o=0 to the completely coherent state with /q—>1 is a second-order-type transition 
under a gradual decrease of the order parameter T/Tiuj. It is not a first-order phase 
transition. A pion gas with r,-ras,T, and N, typical of those in RHIC and LHC, is 
expected to contain a large condensate fraction and a high degree of suppression of 
Bose-Einstein correlation. 

The evaluation of the two-particle correlation function indicates that the usual 
“chaoticity parameter” A can only be interpreted as an experimental tool to label the 
intercept of the correlation function C{p, q) at q—0. The parameter A is correlated 
with the degree of chaoticity only for small values of p but is at variance from such 
an interpretation of chaoticity at high values of p, as shown in Figs. 4 and 5(a). 

We have written out the functional form of the three-particle distribution func¬ 
tion as a function of the momenta of the three particles that contains the proper 
chaotic and coherent limits. It permits the description for the transition from the 
chaotic states to coherent states. These results will allow the evaluation of the 
three-particle correlation function in an exactly solvable problem that will assist the 
comparison with three-particle correlation measurements. 
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